Charge Pumping in Mesoscopic Systems coupled to a Superconducting Lead 
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We derive a general scattering-matrix formula for the pumped current through a mesoscopic 
region attached to a normal and a superconducting lead. As applications of this result we calculate 
the current pumped through (i) a pump in a wire, (ii) a quantum dot in the Coulomb blockade 
regime, and (iii) a ballistic double-barrier junction, all coupled to a superconducting lead. Andreev 
reflection is shown to enhance the pumped current by up to a factor of 4 in case of equal coupling 
to the leads. We find that this enhancement can still be further increased for slightly asymmetric 
coupling. 
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I. INTRODUCTION 

It is well-known that the quantum transport proper- 
ties of a mesoscopic system are modified in the presence 
of a superconducting interface, due to interference be- 
tween normal and Andreev reflections. Andreev reflec- 
tion (AR)El is the electron-to-hole reflection process which 
occurs when an electron with energy slightly above the 
Fermi energy is incident on the boundary between a nor- 
mal metal and a superconductor: the electron enters the 
superconductor after forming a Cooper pair, and leaves a 
hole in the normal metal with energy slightly below the 
Fermi level which travels back along (nearly) the same 
path where the electron came from. Because of the phase- 
coherent character of AR, it is interesting to study its ef- 
fect on transport in mesoscopic systems, where phase co- 
herence plays an important role. In the last decade, this 
has led to the discovery of a wealth of quantum interfer- 
ence effects in mesoscopic normal-metal-superconductor 
(NS) structures^, such as the observation of a large nar- 
row peak in the differential conductance of a disordered 
NS junction ("reflectionless tunneling" )□, and the discov- 
ery of novel Kondo phenomena in quantum dots attached 
to a normal and a superconducting leado. In addition, 
investigations of the conductance in superconductor- 
carbon-nanotube devices have recently appearedu, which 
indicate that also in these devices resonant behavior due 
to AR occurs. 

The purpose of this paper is to study the effects of AR 
on a new type of mesoscopic transport, namely adiabatic 
quantum pumping. Quantum pumping involves the gen- 
eration of a d.c. current in the absence of a bias voltage 
by periodic modulations of two or more system parame- 
ters, such as e.g. the shape of the system or a magnetic 
field. The idea was pioneered by Thouless for electrons 
moving in an infinite one-dimensional periodic potentials. 
In recent years, adiabatic quantum pumpjjia-. in quan- 
tum dots has attracted a lot of attentionQulj. Quan- 
tum dots are small metallic or semiconducting islands, 
confined by gates and connected to electron reservoirs 
(leads) through quantum point contacts (QPCs)£3. In 
addition to investigations of pumping in quantum dots, 
theoretical ideas have been put forward for charge pump- 



ing iB-,carbon nanotube and for pumping of Cooper 
pairaia. Here we consider a mesoscopic system consisting 
of an arbitrary normal-metal region, e.g. a quantum dot 
or QPC, coupled to a superconductor, as schematically 
depicted in Fig. |l|(a) . We start by deriving a general for- 
mula for the pumped current through this NS system in 
terms of its scattering matrix (Sec. II). This is the N- 
mode generalization of the result of Wang et al. for a NS 
system with single-mode leadstU. We then use this result 
to calculate the c urrent pumped through: a simple peri- 
staltic pump (Sec. niA ), a quantum dot in the Coulomb 
blockade regim e (Sec. [II B ) , and a double-barrier junc- 
tion (Sec. Ill Cj ), each coupled to a superconducting lead. 
Comparing these with the pumped current in the corre- 
sponding systems attached to normal leads only, shows 
that AR enhances quantum pumping by up to a factor of 
> 4 in systems with (nearly) symmetric coupling to the 
leads, while it reduces quantum pumping in the opposite 
situation of strongly asymmetric coupling. 



II. DERIVATION OF NS PUMPING FORMULA 

Consider the system in Fig. |l](a). The normal region 
(which may e.g. be disordered, or contain a constric- 
tion) is coupled to ideal normal leads 1 and 2 containing 
N modes each. No bias voltage is applied to the sys- 
tem, so all reservoirs are held at the same potential. We 
assume a constant pair potential A(r) = Aoe 1 ^ in the 
superconductor, which is applicable for wide junctionstfl 
and has previously beep-used to derive the conductance 
through a NS junctionEJ. We also assume that the NS 
interface is ideal, i.e. no specular reflection occurs for 
energies < e < Aq, with e the energy measured from 
the Fermi energy e^. The. scattering matrix 5ns of the 
entire system is given byE3 



Sns(c) 



S ee {e) S eh (e) 
S he {e) S hh (e) 



(1) 



where S ee -S hh are N x N scattering matrices given by 



S* ee (e) 
S eh (e) 
S he (e) 
S hh (e) 



r u (e) + a 2 t l2 {e)r* 22 {-e)M e t 2l {e), (2a) 
ae^t 12 (e)M h t* 21 (-e), (2b) 
ae- 4 **J a (-e)M e t 21 (e), (2c) 
r^-e) + a 2 ^ 2 (-e)r 22 (e)A4^ 1 (-e).(2d) 



Here r ri (e) [ty(e)], i,j = 1,2, denotes the reflec- 
tion [transmission] amplitude for electrons at energy e, 
with < e < Ao, traveling from lead i [j] to lead i, 
a ee exp[-zarccos(e/A )], M e ee [1 - a 2 r 2 2(e)r2 2 (-e)] -1 , 
and Mh ee [1 — a 2 r 22 (— e)r22(e)] _1 - The scattering ma- 
trix ([I]) is unitary and satisfies the symmetry relation 
Sns ( e , B, = S-Ns(e->—B,—(f))ji for time reversal in- 
variance. Adiabatic quantum pumping in this NS junc- 
tion is obtained by slow and periodic variations of two ex- 
ternal parameters X\ and X 2 as X\ (t) — Xi+SXx sin(wi) 
and X 2 (t) = X 2 + 5X 2 s'm(ujt + <fi). The frequency u> has 
to be such that lu <C TjJ on , with T<j we ii the time parti- 
cles spend in the system, in order for equilibrium to be 
maintained throughout the entire pumping cycle. The 
net charge SQ(t) emitted into lead 1 due to the modu- 
lations 8X\ and SX 2 consists of the amount of negative 
charge carriers (electrons) minus the amount of positive 
charge carriers (holes) emitted into lead 1. For infinites- 
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FIG. 1: (a) Normal-metal region (hatched) adjacent to a su- 
perconducting lead (S). The normal leads 1 and 2 contain N 
modes each, (b) NS pump in a ID wire containing a tunnel 
barrier 7 and an external potential U. (c) Quantum dot cou- 
pled via tunneling barriers Ti and T2 to a normal (1) and a 
superconducting (S) lead. 
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1,2, this charge is given by 
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where the indices a and /? are summed over all N modes 
in lead 1. This expression is obtained along the same 
lines as the pumpedj-current in a quantum dot coupled 
to two normal leadaD and based aa a formula derived 
by Biittiker et a/ED, see also Ref.E3. The total charge 
emitted into lead 1 during one period r ee 2tt/lu is found 
by integrating Eq. (||) over time, 
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and rewriting this as an integral over the area A that 
is enclosed in parameter space (Xi,X 2 ) during one pe- 
riod. We then find that the total current Ins = <f~*2( T ) 
pumped into lead 1 is given by 



h 



NS 



J ^^dx x dx 2 n aP (x 1 ,x 2 



uie 
2 



a,/3£l 



— <5AW 2 sin0 V U aP (X l7 X 2 ), 

ZTT * — ' 



(4a) 
(4b) 



with 



n Q/ 3 (Xi , X 2 



Im 



dX 1 dX 2 



fi Qhe * Qhe 

dX 1 dX 2 



(5) 



Eq. (Q) is valid at zero temperature and to first order 
in the frequency cui— For N—l it reduces to the single- 
mode result of Ref.li3 (apart from a factor of 2 for spin 
degeneracy in the latter). Eq. ( |4b| ) applies for bilinear 
response in the parameters X% and X 2l in which case 
the integral in (4a) becomes independent of the pumping 
contour. Jjjs is of similar generality as the expression for 
the piumped current in the presence of two normal-metal 
leadsQ 



E 



Im 



aGl /3£{1,2} 



9S a/3 dS a fj 
dX 1 dX 2 



• (6) 



Here S a p denotes the 2N x 2N scattering matrix of the 
system. Note that the index /3 in this case is summed 
over the modes in both lead 1 and lead 2, since electrons 
can be incident from either lead. In our NS junction, 
the charge pumped into the right lead is converted into 
a supercurrent. 

In order to illustrate the result (^), we now proceed to 
apply it to several NS configurations. Unless otherwise 
noted, we restrict ourselves to the linear response regime 
corresponding to weak pumping. In that regime, only 
the scattering matrix at the Fermi level e = is needed. 
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III. APPLICATIONS 
A. Peristaltic pump in a one-dimensional (ID) wire 

As a first example, consider a "peristaltic" pump 
formed by a ID wire in contact with a superconduct- 
ing lead, see Fig. |l|(b). The pump is operated by peri- 
odically opening and closing a tunnel barrier of height 
7 and varying an external potential U. Calculating the 
scattering matrix (|l|) from the Schrodinger equation with 
potential Vjx) — 7? (x) + UQ[x(L — x)] and substituting 
into Eq. (f4b|) yields the current 



3V2 



'NS 



32 7rfc ; 



5U sin 2 (fcL). 



(7) 



Comparing this with-the analogous expression in the case 
of two normal leadsLJ, 



ujeL 
8ir 2 k 



SU 



167T 2 fc 2 



6U(irsm 2 (kL) - sin(2fcZ)), (8) 



we see that the presence of AR can both enhance and 
reduce the pumped current. Maximum enhancement of 
Ins /In — 3v2/2 ~ 2.1 occurs for short wires such that 
2k L <C 1. Note that 7ns > unlike In, has no classical 
contribution and is entirely due to quantum interference 
in the wire. 



B. Pumped current peak heights in a 
nearly-isolated quantum dot 

As a second example, we study quantum pumping 
in a quantum dot which is coupled via two tunneling 
barriers with transmission probabilities T± and T 2 to 
a normal and a superconducting single-mode lead, sec 
Fig. 0(c). Pumping is achieved by periodic variations of 
the strength of the two tunneling barriers V± and V 2 (with 
a/(1 — Ti :2 )/Ti i2 for delta-function barriers) as 



v> 



1,2 = 

V x {t) = Vx+SVi sm(ut) and V 2 (t) = V 2 +5V 2 sin(wi- 
We are interested in the regime of high barriers (where 
transmission is low, T\,T 2 <C 1), and weak pumping 
(8V m <§; V m ,m — 1,2). At low temperatures such that 
ksT « A< Ec [with A the single-particle level spacing 
and Ec=e 2 /C the charging energy of the dot, G being 
the total capacitance] the quantum dot then remains in 
the Coulomb blockade regime during the whole pump- 
ing cycle and transport through the dot is-tnediated by 
resonant transmission through a single^jtiellMj. Substitut- 
ing the appropriate scattering matrixr 9 rr°l into Eq. ( |4b| ) 
yields, up to lowest order in T± and T 2 and for thermal 
energies less than the total decay width T into the leads 
k B T < T < A [with T = Ti + T 2 = Hv(Ti + T 2 ) and 
v the attempt frequency, the inverse of the round-trip 
travel time between the two barriers! , 



uje 

Ins = -7— 



^-8V\5V 2 sin< 

47T 



,7/2, 



[\(T 2 + T 2 ) 



/ 7Y + V7¥) 



(9) 



Here e rcs denotes the resonance energy for a completely 
isolated dot (T x =T 2 =0). Note that Eq. (|) is not sym- 
metric with respect to T\ and T 2 , in contrast with 



the conductance^ Gns = it 



3spe 

60! 



T?T? 



- [iP ? +I !) + (^)"]. through 
this system. This is due to the fact that Ins depends on 
dS/dX, whereas Gns depends on the transmission eigen- 
values, the eigenvalues of the matrix t\ 2 ti 2 . Compared 
to th 
leads 



.pumped current in a dot coupled to two normal 



In = —6V\6V2sm4> : : — ; -5 (J^O 

\\{Ti + T 2 y + (^) 2 Y 



we find that 

Ins _ 
1n~ ~ ~ 



T 2 T 2 



\(T 1+ T 2 f 



(11) 



(Ti=T 2 ) close to resonance, 
AR enhances the pumped 



For symmetric coupling 
(e - e rcs ) < hvT lt2 
current by a factor of 4EH. In case of strongly asymmet- 
ric coupling close to resonance, on the other hand, AR 
reduces the pumped amplitude: Ins /In ~ {Ti/T 2 ) 2< ^ 
for Ti<T 2 , and Ins/In - T 2 /Ti<l for T 2 «.Ti. The 
enhancement by a factor of 4 for symmetric barriers con- 
sists of 2 contributions of a factor of 2: one factor of 2 
is due to the contribution of both electrons and holes to 
the current, which is also responsible for the doubling of 
conductance Gns/Gn=2 in this NS structuralj. A sec- 
ond factor of 2 comes from the asymmetry of the NS 
dot with respect to injection of charge carriers into the 
leads, since electrons can only leave the system through 
the left, normal lead. This leads to an extra doubling of 
the pumped current compared to the normal case where 
electrons are injected into both the left and the right 
leads. This extra factor of 2 does not occur in the pres- 
ence of an applied bias, as for conductance, since the bias 
causes charge carriers to flow from one side to the other 
in both the normal and the NS system. Note that due to 
the asymmetry of Eq. ([ll]) with respect to T% and T 2 the 
maximum attainable enhancement is even larger than 4: 
for a slightly asymmetric junction (with Ti/T 2 ^1.26) one 
obtains Ins/In^ 4.23. In this case quantum interference 
between electrons and holes in the NS system is maximal. 
If the barrier asymmetry is further increased, Jns /In de- 
creases and eventually becomes less than 1 for strongly 
asymmetric coupling, when pumping is dominated by one 
barrier only. 

At temperatures higher than the decay width, 
r^fc^r^A, the pumped current exhibits Coulomb 
oscillations as a function of an applied gate 
voltaga!3. The peak heights of these oscillations 
can be obtained by thermally averaging Eq. (^) as 



^NS,pcak = — / delftsf (t,T) 
where f(e,T)=[l + cxp(e/fc B T)] 



ik B T ■ 



denotes the Fermi 



4 



function. We obtain 



<NS,pcak 



167rfcsT 



(I? 



7?) 



This thermal average does not explicitly include the ef- 
fect of the charging energy Eq on the pumped current. 
A full linear response theory for Coulomb blockade con- 
ductance oscillations including charging energy was de- 
veloped in RefJla. There it was shown that for temper- 
atures ksT <C A only one level iV m i n participates in the 
transport [iV m i n is defined as the level which minimizes 
the energy J5n+U(N)-U(N-1)- ep, with E^ the energy of 
the Nth level of the dot, and U(N) the electrostatic en- 
ergy of a dot containing N electrons] , and the oscillation 
peaks are well described by the thermal average. One 
can show that for the same reason the pumped current 
peaks in this temperature range are well described by the 
thermal average (|l2|), with the understanding that T\ <2 
in Eq. (|l2] ) refer to the level JV m i n . From ( |l2] ) and the 
analogous normal-state resultta 

iN.poak = 1g _,„ SV 1 dV 2 sm0 (rr \ ^ ^ 2 (13) 



m + r 2 )2 



we obtain 



/NS.pcak _ /8 + 3tt\ TlT 2 {T x + T 2 



'N.pcak 



• T|)5/2 



(14) 



Also here, AR enhances the pumped current in case 
of symmetric tunnel barriers, while a reduction occurs 
for asymmetric barriers. Maximum enhancement of 
^NS,pcak/^N, P oak ~ 2.55 is reached for Ti/T 2 ~ 1.292. This 
factor is less than 4, because the average over energy e 
from which the peak heights (|lj) are obtained also in- 
volves contributions of Ins [Eq. (^)] further away from 
resonance, for which /ns/7m is much less than 4 [consider 
e.g. Eq. (0) for e - e res = \%v yj T% + T%\. This results 
in lower maximal enhancement of the pumped current 
peaks (|l4|) at higher temperatures fcflT> V. 

Another interesting result is obtained in this system 
by relaxing the assumption of weak pumping and con- 
sidering quantum pumping by varying the two tunnel- 
ing barriers in such a way that the loop which describes 
the pumping cycle in parameter space encircles the en- 
tire resonance line. For normakmctal contacts this prob- 
lem has recently been studiedEJ and led to the predic- 
tion that at zero temperature the charge transferred dur- 
ing one pumping cycle is quantized, Q = e (for spin- 
less electrons). The transferred charge in our NS system 
is obtained by substituting the scattering matrix of a 
ID double-barrier junction given in RefJij into Eq. (^a) 
and integrating over the resonance line V^ 1 + V 2 X = 
|(e-e rcs )/H « 1, with Vr 1 = ^/T,/(l - T t ) for % = 1, 2. 
Wc then obtain 



3V2e 



/ Z [l 



(l + z) 2 (l-z 2 ) 



2\:>, 



Q Z 2 + z 4]5/2 



(15) 



so AR neither enhances nor reduces quantum pumping. 
This occurs because charge is effectively transferred by 
a shuttle mechanism (first through one barrier and then 
through the next), which is unaffected by Andreev in- 
terference effects and fixed by the pumping loop. Since 
only pairs of electrons can enter the superconductor, but 
the strong Coulomb interaction (charging energy) forbids 
simultaneous pumping of 2 electrons with opposite spin, 
this pumping process is not allowed in a_nearly-closed 
NS quantum dot. As pointed out in Ref.Ej it can, how- 
ever, occur in a double-barrier junction in which electron- 
electron interactions may be neglected. Both in case of 
two normal and in case of one normal and one supercon- 
ducting single-mode lead a quantized amount of charge 
of 2e is then transferred during each pumping cycleEa. 



C. Ballistic double-barrier junction 

Finally, we compare the pumped current In i n a 
NIiNI^N junction vs. Ins in a NI1NI2S junction, 
where Ii <2 denote tunnel barriers with transmission prob- 
ability per mode Ti 2 , see inset of Fig. p[ In linear re- 
sponse, assuming iV-mode leads and ballistic transport 
between the barriers, the pumped currents are given by 



7 N = IcT, T 2 ^ — 



and 



Ins = 4IcT?Ti 



n=l 



N A 

n=l 



„ , C N s'mcj) n 
(£>n + £ncos<M 2 ' 



j- B NS cos <f) n + C NS sin < 
{Dns + Ens cos0 n ) 3 



■-(17) 



2tt 



6V16V2 sin 



and A 



N. 



ns are given 



V R 2T1T2 



Here Iq 

by A N = jTT{r 2 + v^TT, B N = -V RjhjVRtfi + 
y / Wl\), C N = -l +R 1 R2-VR 1 R 2 T 1 T 2 , D N = l+R x R 2 , 
E N = A NS = 2^Rl, -Bns = 

VRl(2-T 2 ), C NS = -{JT X {2-T 2 )- 
(1 + Ri)(l + R 2 ), and E NS = For L > X F , 

with Xf the Fermi wavelength, and NTi ^> 1 the current 
is not dominated by a single resonance and the phases 
<j) n are uniformly distributed from to 2tt. Replacing 
the sums in Eqs. ( |l6| ) and ([l?]) by integrals over <f> n then 
yields the average pumped currents 



(In) 



± 1 ± 2 



n(l-R 1 R 2 ) 



[2(-l + i2ii2a- 



y/R 1 R 2 T 1 T 2 ) - n (y/R^ + y/Rtf\) ,(18) 



and 



(7 NS ) = 47 c 



T\T 2 



<Dls 



Ans {E^s 



E 2 m )^ 
■ 2C NS )/2 



2Cns£'ns(7 , ns — ^ns) 2 
h3B NS D NS E NS /2] , (19) 
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FIG. 2: Ratio of the average pumped currents {Ins) and 
(In) in ballistic NININ and NINIS junctions as a function of 
T2 for Xi=l, 0.5, and 0.1. The inset shows the double-barrier 
junction considered. 



whose ratio is plotted in Fig. As for a nearly-isolated 
quantum dot, (Ins) / {In} is largest in case of nearly- 
symmetric coupling. Only for T\ = 1 the maximum en- 
hancement by a factor of > 4 is obtained, since for T\ < 1 
less electrons reach the NS interface due to normal reflec- 
tions at the barriers, which in the absence of resonances 



reduces the effect of Andreev reflection. 



IV. CONCLUSION 

In conclusion, we have studied adiabatic quantum 
pumping in mesoscopic NS systems. Compared to the 
conductance in these systems we predict two striking 
differences: (1) For a nearly-isolated quantum dot with 
symmetric (T±=T2) tunneling barriers, and a transparent 
(Ti=X2=l) double-barrier NS junction, AR enhances the 
pumped current by a factor of 4, which is twice the max- 
imum enhancement of the conductance in these systems. 
(2) In case of quantum pumping this enhancement is not 
an absolute maximum, whereas in case of conductance it 
is. These differences are due to, resp., the absence of an 
external bias and the asymmetric dependence on the tun- 
nel barriers in case of quantum pumping. We hope that 
these fascinating effects of Andreev reflection on quan- 
tum pumping will find experimental confirmation, .-e.g. 
in present-day available nearly-closed quantum dotscj. 
Stimulating discussions with CM. Marcus, N.W.A. Tol- 
lenaar and Yu. V. Nazarov are gratefully acknowledged. 
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menteel Onderzoek der Materie (FOM), by NFS grant 
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